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Abstract The carbon isotopic composition of 
methane produced in anoxic marine sediment is 
controlled by four factors: (1) the pathway of 
methane formation, (2) the isotopic composition of 
the methanogenic precursors, (3) the isotope 
fractionation factors for methane production, and 
(4) the isotope fractionation associated with methane 
oxidation. The importance of each factor was 
evaluated by monitoring stable carbon isotope ratios 
in methane produced by a sediment microcosm. 
Methane did not accumulate during the initial 42-day 
period when sediment contained sulfate, indicating 
little methane production from "noncompetitive" 
substrates. Following sulfate depletion, methane 
accumulation proceeded in three distinct phases. 

First, C0 2 reduction was the dominant methanogenic 
pathway and the isotopic composition of the methane 
produced ranged from -80 to -94 %o. The acetate 
concentration increased during this phase, suggesting 
that acetoclastic methanogenic bacteria were unable to 
keep pace with acetate production. Second, acetate 
fermentation became the dominant methanogenic 
pathway as bacteria responded to elevated acetate 
concentrations. The methane produced during this 
phase was progressively enriched in 13 C, reaching a 
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maximum 5 l3 C value of -42 %o. Third, the acetate 
pool experienced a precipitous decline from >5 mM 
to <20 jiM and methane production was again 
dominated by CO 2 reduction. The 5 13 C of methane 
produced during this final phase ranged from -46 to 
-58 %o. Methane oxidation concurrent with methane 
production was detected throughout the period of 
methane accumulation, at rates equivalent to 1 to 8% 
of the gross methane production rate. Thus methane 
oxidation was too slow to have significantly 
modified the isotopic signature of methane. A 
comparison of microcosm and field data suggests 
that similar microbial interactions may control 
seasonal variability in the isotopic composition of 
methane emitted from undisturbed Cape Lookout 
Bight sediment. 

INTRODUCTION 

Recent estimates place the current global methane 
flux to the atmosphere at 543 ± 95 x 10 12 g CH 4 yr* 1 
[Cicerone and Oremland, 1988]. A relatively high 
level of precision (±17%) is possible because the 
global methane flux is calculated from three well- 
characterized parameters: the size of the atmospheric 
methane reservoir, the rate of change in atmospheric 
methane concentration, and the average residence 
time for atmospheric methane. Accurate estimates of 
reservoir size and rate of change are provided by 
concentration time series data from permanent 
stations throughout the world [Blake and Rowland, 
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1988]. Likewise, the residence time for methane is 
tightly constrained by atmospheric models of the 
hydroxyl radical field (calibrated with methyl 
chloroform distributions) and laboratory measure- 
ments of the rate coefficient for the methane- 
hydroxyl reaction [Prinn et al., 1987; Vaghjiani and 
Ravishankara, 1991; Taylor et al, 1991; Fung et al., 
1991]. 

The contribution of individual sources to the total 
methane flux is much more difficult to constrain. The 
vast number and widespread geographic distribution 
of significant methane sources means that detailed 
studies of emission rates are limited to a relatively 
small number of representative systems. The pre- 
cision of source strength estimates based on direct 
flux measurements is limited by spatial and temporal 
variability and uncertainty in the areal distribution of 
ecosystem types. The uncertainty of the global 
methane flux calculated by summing individual 
source terms is about ±50% [Tyler, 1991]. 

Quantitative knowledge of individual source 
strengths is of critical importance because atmos- 
pheric methane concentrations have more than 
doubled during the past 300 years [Peatman, 1986]. 
Since methane is a radiatively active trace gas, an 
increase in the atmospheric methane reservoir could 
alter the Earth's energy budget and contribute to 
global-scale climate change [Cicerone and Oremland, 
1988]. In addition, methane serves as a sink for the 
hydroxyl radical and therefore plays an important 
role in regulating the oxidizing power of the atmos- 
phere [Cicerone and Oremland, 1988]. Accurate 
source strength estimates are necessary to understand 
the underlying causes of increasing atmospheric 
inventories and to evaluate the efficacy of different 
mitigation strategies. 

The isotopic composition of carbon in the 
methane molecule serves as a natural tracer of 
source. Methane produced during biomass com- 
bustion is enriched in ,3 C relative to that produced by 
methanogenic bacteria [Stevens and Engelkemeir, 
1988]. Fossil methane (from natural gas pipelines, 
coal mining activities, clathrate decomposition, and 
peat remineralization) can be distinguished on the 
basis of its low 14 C content [Wahlen et al., 1989]. 

The isotopic composition of atmospheric methane 
can provide a global-scale estimate of the contri- 
bution of bacterial, biomass burning, and fossil 
sources to the global methane budget [Quay et al., 
1991]. 

The accuracy of source strength estimates based 
on isotopic constraints is limited by variability in the 


isotopic composition of methane produced by 
individual sources. For example, 8 13 C values of 
methane emitted from a particular ecosystem (such as 
a temperate wetland, Arctic tundra, or rice paddy) 
can vary by as much as 20 to 30 %o [Quay et al, 
1988; Tyler et al, 1988; Chanton and Martens, 

1988]. The factors that contribute to this variability 
are not well understood. 

One important factor in controlling the stable 
carbon isotopic composition of bacterial methane is 
the mechanism of methane formation. Acetate fer- 
mentation and CO 2 reduction are thought to be the * 
dominant methane production pathways in nature 
[Oremland, 1988]. However, methanogenic bacteria 
are known to produce methane from alternative 
compounds such as formate, methanol, methylated 
amines, and methylated sulfur compounds 
[Oremland, 1988]. These compounds serve as 
important methane precursors in some environments 
[Strayer and Tiedje, 1978; Oremland et al., 1982; 
King et al., 1983; King, 1988], but their widespread 
significance remains a matter of debate. 

The isotopic composition of methane produced 
by a particular pathway depends on the 8 13 C of the 
methane precursor and the isotope fractionation 
associated with the production process (i.e., the 
kinetic isotope effect). Therefore variations in 
S 13 C-CH 4 are controlled by the relative contribution 
of different methane production pathways, the 
isotopic composition of the methane precursors, and 
the kinetic isotope effect for each pathway. In 
addition, aerobic and anaerobic methane oxidation 
can alter the isotopic composition of the residual 
methane reservoir [King et al., 1989; Alperin et al., 
1988]. 

Most field-based studies designed to characterize 
the isotopic composition of methane emitted from 
specific sources have not conducted the process-level 
research necessary to understand the factors that 
control 8 13 C-CH 4 variability. In this study, we moni- 
tored the 8 13 C of methane produced by a laboratory 
microcosm containing anoxic marine sediment from 
Cape Lookout Bight, North Carolina. We concur- 
rently measured methane production rates from C0 2 
and acetate, methane oxidation rates, 8 I3 C-ZC0 2 , 
and a suite of key indicators of microbial processes. 
We observed large variations over time in the 8 13 C of 
methane produced by the microcosm. These vari- 
ations can be attributed both to shifts in the relative 
importance of C0 2 reduction and acetate fermen- 
tation, as well as changes in the isotopic composition 
of the methane precursors. The variability in 8 l3 C- 
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CH 4 can be understood in terms of syntrophic and 
competitive interactions between the microorganisms 
involved in organic matter remineralization. A 
comparison of microcosm and field data suggests 
that similar microbial interactions may control 
seasonal variability in the isotopic composition of 
methane emitted from Cape Lookout Bight sediment. 


METHODS 

Study Site and Sediment Sampling 

Sediment was collected from Cape Lookout 
Bight, a shallow (10 m) barrier island lagoon located 
70 km southwest of Cape Hatteras, North Carolina. 
Organic carbon deposition at this site (165 ± 20 mol 
m‘ 2 yr' 1 ) is among the highest reported for a coastal 
marine environment [Martens et al., 1992]. Sediment 
oxygen penetration is restricted to several hundred 
microns despite an air-saturated water column mixed 
by strong tidal currents. Sediment metabolism is 
dominated by microbial sulfate reduction and 
methane production, with methanogenesis 
accounting for about 30% of the annual integrated 
organic matter remineralization [Martens and Klump, 
1984], Metabolic rates exhibit strong seasonality 
driven by annual temperature changes of >20°C 
[Klump and Martens, 1989]. Sulfate concentration 
profiles are modulated by seasonal variations in 
sulfate reduction rates; the vertical extent of the 
sulfate reduction zone oscillates from >30 cm in the 
winter to <10 cm in the summer [Crill and Martens, 
1987]. Methane production rates are highest 
immediately below the sulfate reduction zone [Crill 
and Martens, 1986] and sediments are methane- 
saturated below this horizon. During the summer, 
methane bubbles form below 10 cm and low tide 
ebullition occurs in response to reduced hydrostatic 
pressure [Martens and Klump, 1980]. The 5 13 C of 
methane contained in these bubbles exhibits seasonal 
variations of ~10 %o [Martens et al., 1986; Burke et 
al., 1987]. Sediment accumulation rates at this site 
are approximately 10 cm yr' 1 [Chanton et al., 1983; 
Canuel et al., 1990]. 

Sediment samples used in this study were col- 
lected in July 1990 when sulfate reduction and 
methane production rates were at their annual max- 
imum. Cores were collected by SCUBA diver in 
22 cm x 13 cm diameter Plexiglass tubes. The tubes 
were gently inserted into the sediment and capped on 
top and bottom before being carefully withdrawn. 


The cores were kept in the dark at in situ temperature 
(26.5°C) during transport to the laboratory. 

Incubation Vessel 

A large capacity syringe (hereafter referred to as 
the "magnum syringe") capable of long-term 
headspace-free incubations was constructed for this 
study (Figure 1). The vessel was designed as a 
syringe to allow for intermittent sample withdrawal 
without introducing a headspace. The presence of a 
headspace would significantly alter concentrations of 
dissolved methane and ZCO 2 . 

The syringe barrel was fabricated from 4-inch 
(10.16 cm) ID Kimax Beaded Process Pipe with a 
10-mm bore stopcock (Teflon plug) welded to one 
end. The piston was constructed of solid polyvinyl 
chloride (PVC) machined to provide a slide-fit within 
the syringe barrel. A thin Teflon plate was secured to 
the bottom of the piston to prevent contact between 
sediment and PVC. Attempts at making a direct 
O-ring seal between the piston and the barrel failed 
due to variability in the inner diameter of the glass 
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Fig. 1. Magnum syringe for incubating ~2 L 
sediment under headspace- free, anoxic conditions. 
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pipe. Therefore the piston was sealed to the barrel by 
means of a custom fabricated coupling that provided 
an O-ring seal to the barrel flange and a double 
O-ring seal around the piston. Sediment within the 
magnum syringe could be homogenized without air 
contact by churning with a perforated Teflon mixing 
plate connected to a Teflon rod that passed through a 
slightly oversized hole bored through the center of 
the piston. The mixing rod was sealed to the top of 
the piston with an Ultra-torr fitting (Cajon Vacuum 
Products) that provided a gas-tight seal while 
allowing for axial movement. 

The magnum syringe was filled with sediment 
from the 0 to 3 cm depth interval. Six cores were 
pooled to give a total sediment volume of ~2 L. The 
piston was positioned in such a way as to exclude 
gas pockets and fastened in place by securing the 
coupling. The magnum syringe was sampled at 
weekly intervals by displacing 120 mL of sediment 
into two 60-mL catheter-tip plastic syringes that fit 
snugly in the arm of the stopcock. A large chain 
clamp immobilized the piston between sampling 
intervals. The sediment was homogenized prior to 
each sampling by churning sediment through the 
mixing plate. The magnum syringe was incubated in 
the dark at 22 ± 1°C. 

Concentration Measurements 

Sulfate, acetate, and ZC0 2 concentrations were 
determined on interstitial water exttacted from whole 
sediment by centrifugation. Care was taken to pre- 
vent sediment and interstitial water from contacting 
air. Sediment was transferred to an argon-flushed 
30-mL centrifuge tube and spun at 6000 G for 
15 min. The centrifuge tube was filled to the top to 
minimize loss of gaseous C0 2 into the headspace 
volume. Following centrifugation, supernatant from 
just above the solid-liquid interface was pipetted 
using a long stainless steel needle fixed to a 5-mL 
glass syringe. The interstitial water was passed 
through a 0.45-qm syringe filter; the first milliliter of 
filtrate was discarded. 

Sulfate concentrations were determined using a 
Dionex 2010i ion chromatograph. In order to mini- 
mize sulfide oxidation in the sample, 1.0-mL inter- 
stitial water was filtered directly into an argon- 
flushed vial containing 100-|iL 10% ( v /v) HC1 and 
stripped of volatile sulfur compounds by bubbling 
with 0 2 -free argon. The samples were diluted 1:10 
with distilled water and passed thru an ONGUARD- 
Ag pretreatment cartridge (Dionex Corporation) to 


remove chloride. The large quantity of chloride in 
interstitial water (-600 mM) interferes with low-level 
(<1 mM) sulfate determinations; after chloride 
removal the detection limit for sulfate was 5 |iM. 

Acetate concentrations were detenu ined by 
HPLC using precolumn deiivatization with 2-nitro- 
phenylhydrazine [Albert et al, 1992]. The deriva- 
tives were prepared directly in the pore water and 
separated by a C8 reversed phase column with 
buffered, ion-pairing solvent. The detection limit for 
this method is 0.5 qM; the precision is typically 
±5%. 

Interstitial water for XC0 2 analysis was stored in 
3.5-mL Pyrex screw-cap vials filled to just below the 
rim and capped with Teflon-faced silicone septa. 

Loss of C0 2 into the gas phase was prevented by 
displacing the headspace with a 4-mm Pyrex rod 
inserted through a hole in the septum. The ZCQi 
concentrations were determined using the inorganic 
carbon channel of a Shimadzu TOC-5000 total 
carbon analyzer. Samples (25 |iL) were injected 
directly into a gas stripping chamber containing 
25% ( v /v) H3PO4. The CO 2 released upon acidifi- 
cation was quantified by a nondispersive infrared 
detector. 

Whole sediment methane concentrations were 
measured by a headspace equilibration technique 
[Alperin and Reeburgh, 1985]. Whole sediment 
concentrations include methane dissolved in the pore 
water as well as methane bubbles trapped in the 
sediment. Sediment (2.5 mL) was transferred to a 
tared 30-mL serum vial containing 3-mL 0.1 N 
NaOH, sealed with a silicone stopper, and thor- 
oughly vortex-mixed to disperse the sediment and 
equilibrate the gas and aqueous phase. An aliquot of 
the headspace was analyzed for methane by a 
Shimadzu Mini-2 gas chromatograph equipped with 
flame ionization detector. Standards were prepared 
by injecting pure methane into sealed 30-mL serum 
vials containing 3-mL 0.1 N NaOH and 2. 5-mL 
distilled water. Samples and standards were kept at 
the same temperature and had equal volumes of 
aqueous phase and headspace, so corrections for the 
quantity of dissolved methane were not necessary. 
Whole sediment methane concentrations were con- 
verted to pore water units by dividing by sediment 
porosity. 

Rate Measurements 

Sulfate reduction [Albert et al., 1992], methane 
production from acetate [Crill and Martens, 1986], 
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methane production from CO 2 fCrill and Martens, 
1986], and methane oxidation [Alperin and 
Reeburgh, 1985] rates were determined by previ- 
ously described radiotracer techniques. Briefly, 
whole sediment was transferred to gas-tight 2.5-mL 
glass incubation syringes and injected with microliter 
quantities of radiotracer ( 35 S0 4 2 \ [U- l4 C]acetate, 
[ 14 C]-HC0 3 \ or 14 CH 4 ). The samples were incu- 
bated at 22 ± 1°C for a time interval sufficient for 
1 to 10% of the tracer pool to turnover. The incu- 
bation time was generally 24 hours, except samples 
injected with [U- 14 C] acetate which had incubation 
times ranging from 10 min to 27 hours. At two times 
(93 and 107 days), the samples amended with 
[U- l4 C]acetate were incubated sufficiently long that 
40-60% of the tracer turned over. Methane pro- 
duction rates from acetate at these two times are 
considered to be underestimated. 

Sediment samples assayed for sulfate reduction 
rate were "killed" by freezing the incubation 
syringes. For other rate measurements, the incu- 
bation was terminated by extruding the sediment into 
serum vials containing NaOH solution. The serum 
vials were quickly sealed with butyl-rubber or sili- 
cone stoppers and stored frozen and upside-down (to 
retard gas loss) until analysis. Collection and quan- 
tification of radioisotopes was according to published 
procedures (references given above). The rates were 
calculated from the following equation: 
a a 

Rate = [C] srT , 

where [C] is reactant concentration (i.e., sulfate, 
acetate, LC0 2 , or methane), a and A are activities of 
recovered product and added reactant, respectively, 
a is the fractionation factor for the tracer isotope 
( 14 C or 35 S) relative to the dominant natural isotope 
( 12 C or 32 S), and t is the incubation time. The isotope 
fractionation factor was taken to be 1.12 for methane 
production from CO 2 [Blair et al., 1992], 1.06 for 
methane production from acetate [Blair and Carter, 
1992], 1.02 for methane oxidation [Alperin et al., 
1988], and 1.045 for sulfate reduction [J0rgensen, 
1978]. Note that isotope discrimination for I4 C is 
twice that for 13 C [Stem and Vogel, 1971]. 

Stable Isotope Measurements 

The S 13 C-SCC >2 analyses were done on 1.0-mL 
aliquots of filtered interstitial water stored in 3-mL 
serum vials capped with silicone septa. Samples 
were acidified with 0.5-mL 1 M H 3 P0 4 and stripped 
with helium into a gas purification line [Schaff et al.. 


1992]. Air and C0 2 were separated by a GC column 
packed with porous silica beads (Unibeads IS, 
Alltech Associates) and quantified by a thermal 
conductivity detector (Gow Mac Coiporation). The 
C02 peak was collected in a cryogenic trap (-196°C) 
fashioned from six loops of 1/8 inch (0.318 cm) 
stainless steel tubing. The helium was pumped away, 
and the C0 2 was transferred to a Pyrex breakseal 
ampoule for introduction to the isotope ratio mass 
spectrometer. 

For 5 13 C-CH 4 analysis, 10-mL whole sediment 
was transferred to a 60-mL serum vial containing 
5-mL IN NaOH and sealed with a black butyl rubber 
stopper. The sample and base were thoroughly 
slurried and an aliquot of the headspace was trans- 
ferred to a gas-tight syringe and injected into a 
combustion line. The sample was Hushed with 
helium carrier gas through a cryogenic trap (-I96°C) 
to remove C0 2 and H 2 0 and into a combustion tube 
packed with CuO at 780°C. The C0 2 and H 2 0 
resulting from methane combustion were passed thru 
a multiloop trap maintained at -120°C (to remove 
H 2 0), and the C0 2 was collected in a breakseal 
ampoule. 

Carbon isotopic analyses were performed with a 
Delta E ratio mass spectrometer (Finnigan-Mat) 
equipped with a modified low-volume inlet system 
[Hayes et al., 1977]. The precision of the 5 13 C 
analysis was ±0.1 %o. 

RESULTS 

Changes in sulfate, methane, acetate, and XC0 2 
concentrations during the 114-day incubation period 
are summarized in Figure 2. The sulfate concen- 
tration decreased from 18.6 mM at the outset to 
0.19 mM at 42 days (Figure 2a). After 50 days, 
sulfate leveled off at 10 to 20 pM, comparable to the 
threshold concenttation below which sulfate 
reduction stops for pure cultures of Desulfobacter 
postgatei [Ingvorsen et al., 1984]. 

The methane concentration remained relatively 
constant (~0.2 mM) until sulfate was depleted to 
<0.2 mM (Figure 2a). Following sulfate depletion, 
the methane concenU*ation increased, exceeding the 
saturation level (1.2 mM [Yamamota et al., 1976]) 
after 70 days. Since the immobilized magnum 
syringe piston restricted headspace formation, subse- 
quent methane production resulted in a pressure 
buildup. The coupling that sealed the piston to the 
barrel was not capable of containing positive pres- 
sures and methane was inadvertently lost. The 
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Fig. 2. Changes in pore water concentrations during 
the 1 14-day incubation period, (a) Sulfate and 
methane; (b) acetate and ICO 2 . With the exception of 
acetate, error bars represent the standard deviation of 
triplicate analyses. A standard deviation of ± 5% is 
assumed for acetate concentrations. The absence of 
an error bar indicates that the standard deviation is 
smaller than the symbol. The hatched line denotes the 
time that pore water sulfate concentrations were 
depleted to <0.2 mM. 


decrease in methane concentration after 85 days is 
due to gas loss from the magnum syringe rather than 
net methane consumption (see below). 

The acetate concentration ranged from 1 to 9 |iM 
during the initial 42-day period when sediment 
contained appreciable sulfate (Figure 2b). Acetate 
concentrations <10 pM are typical for sulfate- 
containing marine sediment [Ingvorsen et al., 1984] . 
Following sulfate depletion, acetate began to accu- 
mulate in the sediment, reaching a peak concentration 
of 5.3 mM at -80 days. The period of acetate accu- 
mulation was immediately followed by a period of 
dramatic depletion as concentrations fell to 5 jiM in 
just 14 days. 

The XC0 2 concentration increased rapidly during 
the period that sediment contained sulfate (Figure 


2b). Sulfate and ZC0 2 concentrations during the 
initial 42-day period are highly correlated (r 2 = 0.99), 
suggesting that sulfate served as the dominant 
terminal electron acceptor for organic matter reminer- 
alization. The slope of the ICCF versus sulfate plot 
(A[£C0 2 ]/A[S0 4 2 ']) is -1.7 ±0.1, suggesting that 
sediment organic carbon remineralized in the mag- 
num syringe experiment had an average oxidation 
state somewhat less than carbon in pure carbohy- 
drate. 

The rate of increase of IC0 2 tapered off fol- 
lowing sulfate depletion. This is primarily due to a 
delay in ZC0 2 production caused by pooling of car- 
bon in the acetate reservoir. However, the gas loss 
that occurred when methane concentrations exceeded 
saturation may also have contributed to a reduced rate 
of dissolved ZC0 2 accumulation. Gas equilibrated 
with sediment pore waters (pH 7) having a ZC0 2 
concentration of 50 mM will have a C0 2 partial pres- 
sure of 0.2 atm [Stumm and Morgan, 1981]. There- 
fore ZCC^ concentrations in samples collected after 
methane concentrations exceeded saturation should 
be regarded as lower limits. 

Sulfate reduction, methane oxidation, and 
methane production rates are presented in Figure 3. 
Sulfate reduction rates measured by the 35 S tracer 
method are consistent with rates calculated from 
changes in sulfate concentration. Sulfate reduction 
rates fell off rapidly from the outset (Figure 3a), 
suggesting a rapid depletion of the most reactive 
fraction of the organic carbon pool [Westrich and 
Berner, 1984]. Since sulfate concentrations through- 
out most of the sulfate reduction period were well in 
excess of the half-saturation constant for sulfate 
reducing bacteria (0.2 mM [Ingvorsen et al., 1984]), 
it is unlikely that the initial decline in the sulfate 
reduction rate was due to sulfate limitation. 

Methane oxidation rates remained near baseline 
level throughout the period that sediment contained 
sulfate (Figure 3a). Note that methane concentrations 
during the initial 42 days (~0.2 mM, Figure 2a) were 
sufficient to support moderate rates of oxidation. 
Following sulfate depletion, methane oxidation rates 
began to increase, achieving maximum values 
between 80 and 100 days. The methane oxidation 
rate at all times represents a small fraction (1 to 8%) 
of the total methane production rate (see below). 
Methane oxidation rates integrated over the 1 14-day 
experiment account for only 2.4% of the integrated 
methane production. Net methane oxidation was not 
observed during this experiment. 
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Fig. 3. Changes in rates of microbial processes 
during the 1 14-day incubation period, (a) Sulfate 
reduction and methane oxidation; (b) methane 
production from C0 2 and acetate. Rates of methane 
production from acetate at 93 and 107 days are 
considered to be underestimated because 40 to 60% 
of the tracer turned-over during the incubation 
period. The error bars represent the standard 
deviation of rates measured in two to four (generally 
3) replicate samples. The absence of an error bar 
indicates that the standard deviation is smaller than 
the symbol. The hatched line denotes the time that 
pore water sulfate concentrations were depleted to 
<0.2 mM. 


Methane production rates from C0 2 and acetate 
were negligible during the period that sediment con- 
tained sulfate (Figure 3b). The absence of C0 2 re- 
duction and acetate fermentation in the presence of 
sulfate is consistent with the well-established prin- 
ciple that sulfate-reducing bacteria out-compete 
methanogens for H 2 and acetate [Kristjansson et al., 
1982; Schonheit et al., 1982], Following sulfate 
depletion, methane production rates from C0 2 and 
acetate began to increase, but the two methanogenic 
pathways followed very different patterns. Methane 
production from C0 2 increased steadily, reached a 


maximum rate of -100 |iM d‘\ and gradually 
declined. In contrast, methane production from 
acetate began relatively slowly, increased sharply to 

-400 jaM d 1 , and then decreased dramatically to very 
low levels. 

Relatively constant methane concentrations 
(Figure 2a) and very low rates of oxidation in 
sulfate-containing sediment (Figure 3a) suggest little 
methane production from "noncompetitive" sub- 
strates (i.e., compounds not metabolized by sulfate- 
reducing bacteria). Therefore acetate fermentation 
and C0 2 reduction are assumed to be the dominant 
pathways for methane production throughout the 
experiment. Total methane production rates (Figure 
4a) are calculated as the sum of methane production 
from C0 2 and methane production from acetate. 






Fig. 4. Changes in total methane production rates 
and pathways during the 1 14-day incubation period, 
(a) Total methane production rates; (b) fraction of 
methane derived from acetate (F). The error bars 
were calculated by error propagation. The absence of 
an error bar indicates that the standard deviation is 
smaller than the symbol. The hatched line denotes the 
time that pore water sulfate concentrations were 
depleted to <0.2 mM. 
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The fraction of methane derived from acetate (F) 
varied considerably over the course of the experiment 
(Figure 4b). Following the sulfate-methane tran- 
sition, C0 2 reduction accounted for the bulk of the 
methane production. As acetate accumulated in the 
pore waters (Figure 2b), the methanogenic pathway 
shifted and acetate became the major methane pre- 
cursor. Acetate fermentation dominated methane 
production until the acetate reservoir was fully 
depleted. At this point, C0 2 reduction resumed its 
role as the dominant methanogenic pathway. Over 
the course of the experiment, acetate fermentation 
accounted for nearly 60% of the total methane 
production. 

The IC0 2 and methane reservoir experienced 
large changes in isotopic composition over the course 
of the 1 14-day experiment (Figure 5). During the 
period that sediment contained sulfate, S 13 C-IC0 2 
values decreased only slightly (Figure 5a), despite a 
large increase in XC0 2 concentration (Figure 2 b). 




TIME (days) 


Fig. 5. Changes in S 13 C values during the 1 14-day 
incubation period, (a) 8 13 C-XC0 2 ; (b) 8 13 C-CH 4 . 
Samples were not analyzed for 8 ,3 C during the first 
6 days of the experiment. The hatched line denotes 
the time that pore water sulfate concentrations were 
depleted to <0.2 mM. 


This suggests that the C0 2 produced by organic 
matter remineralization had a S 13 C similar to the 
initial value of the ZC0 2 pool. Following sulfate 
depletion, the ZC0 2 reservoir became progressively 
enriched in 13 C. 

The 8 13 C-CH 4 values remained nearly constant 
(-59 to -61 7 oo) prior to the onset of methanogenesis 
(Figure 5b). The low methane concentrations during 
this initial period reflect methane that was present in 
the 0 to 3 cm interval at the time the sediment was 
collected. The initial 8 , 3 C-CH 4 values are similar to 
those for methane in gas bubbles collected from Cape 
Lookout Bight sediment in July (-59 to -62 7oo) 
[Martens et al., 1986; Burke et al., 1987]. Following 
sulfate depletion, the methane reservoir underwent 
substantial variations in S 13 C. The S 13 C-CH 4 initially 
dropped from -61 to -88 7oo, gradually increased to 
-47 7oo, and decreased slightly toward the end of the 
experiment. 

DISCUSSION 

The isotopic composition of methane in the 
magnum syringe varied by more than 40 7oo over the 
course of the experiment (Figure 5b). In the fol- 
lowing section, we evaluate the relative importance 
of changes in methanogenic pathway, precursor 
isotopic composition, and kinetic isotope effects in 
controlling the observed variations in 8 13 C-CH 4 . 
Since integrated rates of methane production are 
-40 times greater than methane oxidation, it is 
unlikely that methane oxidation played an important 
role in controlling the isotopic composition of 
methane. 

The discussion is divided into four sections. The 
first section examines the relationship between 
8 l 3 C-CH 4 and the methane production pathway. The 
second section provides constraints on variations in 
the isotopic composition of the major methane pre- 
cursors (CO 2 and acetate). The third section de- 
scribes the microbial interactions that may have 
played a role in regulating both methane production 
pathways and substrate isotopic composition. The 
final section compares the results of the microcosm 
experiment with in situ seasonal variations in isotopic 
composition of methane emitted from Cape Lookout 
Bight sediments. 

8 13 C-CH 4 Versus Methane Production Pathway 

Measured 8 13 C-CH 4 values (Figure 5b) reflect 
the isotopic composition of methane accumulated in 
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the magnum syringe. Before we can examine the 
relationship between 5 13 C-CH 4 and methane pro- 
duction pathway, the 5 I3 C of the methane produced 
during a discrete time interval must be determined by 
an isotope mass balance calculation. Since the pos- 
itive pressures that developed in the magnum syringe 
when methane concentrations exceeded saturation 
resulted in gas loss, methane concentrations must be 
corrected before mass balance constraints can be 
applied. 

Methane concentrations corrected for gas loss . 
Methane concentrations calculated by integrating 
production rate data are unaffected by gas loss from 
the magnum syringe and therefore reflect the total 
quantity of methane produced during the experi- 
mental period. The methane concentration at time t 2 is 
equal to the concentrauon at bme ti plus the time 
integral of the total methane production rate (MPR) 
(Figure 4a): 

2 

[CH 4 ] 2 = [CH 4 ], + J MPR dt . (1) 

1 

Calculated (equadon (1)) and measured concen- 
tradons began to diverge as the methane concen- 
tradon approached saturadon (Figure 6). By the end 
of the experiment (1 14 days), the measured methane 
concentration represents less than 10% of the 
calculated value. 

Isotopic composition of methane produced. The 
isotopic composition of methane produced in the 
magnum syringe (8 13 C-MP) during any discrete time 
interval (e.g., t\ to t 2 ) can be calculated by an isotope 
mass balance equation: 

5 13 C-MP = 

(5 13 C-CH 4 ) 2 [CH 4 ] 2 - (S 13 C-CH 4 ), [CH 4 ], _ 

~ [CH 4 ] prod ' W 

where [CH 4 ] PRO d is the concentration of methane 
produced during the dme interval q to t 2 . 

Two factors must be considered before applying 
(2) to the magnum syringe experiment. First, the gas 
phase that appeared when methane concentrations 
exceeded saturation could lead to isotopic fraction- 
ation if there was an isotope effect associated with 
transfer of dissolved methane to the gas phase. 
Bernard et al. [1976] have shown that there is little 
isotopic fractionation between CH 4 (aq) and CH 4 (g) 



Fig. 6. Changes in measured and calculated methane 
concentrations during the 1 14-day incubation period. 
The method used to calculate methane concentrations 
is described in the text. The arrow on the concentra- 
tion axis marks the dissolved methane concentration 
equivalent to a partial pressure of 1 atm. The hatched 
line denotes the time that pore water sulfate concen- 
trations were depleted to <0.2 mM. 

(0.3 %o), so formation of a gas phase would not 
have an appreciable effect on the isotopic compo- 
sition of the dissolved methane pool. 

Application of (2) is further complicated by the 
large fraction of the methane production that was lost 
from the magnum syringe when positive pressures 
developed (Figure 6). Since effusive transport dis- 
criminates between 12 C- and 13 C-methane according 
to Graham's law [Mason and Kronstadt, 1967], gas 
loss via this mechanism would cause methane re- 
maining in the magnum syringe to become enriched 
in 13 C. In order for effusive transport to occur, 
methane must pass through an orifice without col- 
liding with other gas molecules. This requires that 
the diameter of the leak in the O-ring seal be less that 
the mean free path of the methane molecule (46 nm). 
The effusive flux through a hole this size can be cal- 
culated [Mason and Evans, 1969] by assuming the 
leak to have cylindrical geometry with length equal to 
the O-ring cross section (3 mm). The effusive flux 
(~10' 8 mmol d' 1 ) is 7 orders-of-magnitude less than 
the actual leak rate (0.1 to 0.4 mmol d' 1 ) estimated 
from the difference between calculated and measured 
methane concentrations (Figure 6). This implies that 
pressure-driven gas transport from the magnum 
syringe must have occurred by other mechanisms. 
Since there is no mass discrimination during turbu- 
lent or viscous flow [Halsted and Nier, 1950], gas 
loss through faults in the O-ring coupling is unlikely 
to alter the 5 l3 C value of the methane. 
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Methane loss did reduce the size of the reservoir, 
however, so that the isotopic composition ot any 
methane produced had a proportionately greater 
effect on the cumulative 5 13 C. The effect of gas loss 
on the isotopic composition of the methane reservoir 
was corrected in the following way. Cubic spline 
functions [Ahlberg et al., 1967] were fit to the meas- 
ured and calculated methane concentration (Figure 6) 
and 5 13 C-CH 4 data (Figure 5b) to allow interpolation 
at any q and t 2 . Loss of methane from the magnum 
syringe was compensated by calculating the methane 
concentration at t 2 as 

[CH 4 ] 2 = [CH 4 ]i + [CH 4 ] PROD » (3) 

where [CH 4 ]prod = [CH 4 ] 2 - [CH 4 ] t . (4) 

The asterisks indicate calculated (as opposed to 
measured) methane concentrations. Thus [CH 4 ] 2 in 
(3) represents the methane concentration that would 
be present at t 2 if there was no gas loss during the 
time interval ti to t 2 . Any gas loss is treated as though 
it occurred during an infinitesimal time interval just 
after t 2 . 

The 8 l3 C of methane produced during the time 
interval q to t 2 was corrected for the effect of gas loss 
by combining (2), (3), and (4). The time step (t 2 - q) 
was made sufficiently small (0. 1 day) that the calcu- 
lated 8 13 C-MP closely approximates the 8 13 C of 
methane produced during any instant. The 8 13 C value 
of methane produced during a finite time interval 
(e.g., between consecutive magnum syringe 
samples) was then calculated by integrating 
instantaneous [CH 4 ] P rod and S 13 C-MP values: 

B 

J([CH 4 ] P rod5 i 3 C-MP) dt 

8 l3 C-MP = p- , (5) 

|[CH 4 ]prod dt 

A 

where A and B represent the beginning and end of 
the time interval of interest. 

The 8 l3 C of methane produced by the sediment 
microcosm (calculated by equation (5)) varied 
between -94 and -42 °/oo (Figure 7). The range 
(52 %o) is greater than that of the measured values 
because the cumulative reservoir represents an inte- 
grated isotopic composition. The isotopic compo- 
sition of the cumulative reservoir tracks 
8 l3 C-MP fairly closely (cf. Figures 5b and 7) 
because methane continuously escaped from the 


magnum syringe once the methane concentration 
exceeded saturation. 

There is no simple relationship between the 8 l3 C 
of the methane produced in the magnum syringe and 
the methane production pathway. For example, C0 2 
reduction dominated methane production (accounting 
for -70% of total production) both at the outset of the 
methanogenic period and toward the end of the ex- 
periment (Figure 4b). However, methane produced 
during the initial period (-80 to -94 °/oo, Figure 7) 
was highly depleted in 13 C relative to the methane 
produced during the later period (-46 to -60 %o). A 
simple relationship between 8 13 C-CH 4 and methane 
production pathway is expected only if the methane 
produced by each mechanism has a distinct isotopic 
composition. The absence of such a relationship 
indicates variability in the 8 l3 C value of methane 
produced by a particular pathway. Such variability 
can be driven by changes in the isotopic composition 
of the methane precursors (C0 2 and acetate) and/or 
changes in the magnitude of the kinetic isotopic effect 
associated with each methane production mechanism. 

Overall isotope mass balance. Conservation of 
mass requires that the end products of organic matter 
degradation (C0 2 and methane) have the same cumu- 
lative isotopic composition as the remineralized 
organic carbon. The isotopic composition of the 



Fig. 7. Changes in isotopic composition of methane 
produced (S 13 C-MP) during the 1 14-day incubation 
period. The solid horizontal bars represent the S 13 C 
of methane produced between consecutive sampling 
intervals. The shaded region represents the 8 l3 C of 
methane produced by C0 2 reduction. Details of the 
calculations are given in the text. The hatched line 
denotes the time that pore water sulfate concen- 
trations were depleted to <0.2 mM. 
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labile organic carbon should not be directly affected 
by the nature of the terminal electron acceptor. There- 
fore the 8 I3 C of organic carbon remineralized during 
the period of sulfate reduction (8 ,3 C-C 0 r G [SR]) 
should be similar to the 5 13 C of organic carbon 
remineralized during the period of methane pro- 
duction (8 i3 C-C 0 rg[MP]). 

Since CO 2 is the sole end product of organic 
matter decomposition in sulfate-containing sediments 
(Figure 2), the 5 13 C of organic carbon remineralized 
during the period of sulfate reduction is simply calcu- 
lated as the 8 13 C value of the XCO 2 produced: 

5 i3 c c [SR] = (8 13 C-ZCQ 2 ) f1n [£CQ 2 1p, m 

ORG [SC0 2 ]p [N - [IC0 2 ] INI 

. (8 13 C-ZC0 2 ), m [ZC 0 2 ], n , 

[£C0 2 ] F{N - [SCO 2 ] INI 

where INI and FIN denote [XC0 2 ] or 8 13 C-IC0 2 at 
the start of the experiment and at the time that sulfate 
became depleted (day 42), respectively. Between 14 
and 42 days, the ZC0 2 concentration increased from 
37.6 to 55.2 mM (Figure 2b), while the S 13 C-IC0 2 
decreased from -13.4 to -14.1 %x> (Figure 5a). The" 
8 13 C value of organic carbon remineralized during 
the period of sulfate reduction (S 13 C-C 0 r G [SR]) is 
-15.6 %o (equation (6)). This value differs from the 
S 13 C of the bulk organic carbon in Cape Lookout 
Bight sediments (-18 to -20 %o) [Blair et al, 1987; 
Martens et al., 1992; Blair and Carter, 1992], 
suggesting preferential remineralization of an organic 
matter component that is enriched in l3 C relative to 
the bulk phase. This result seems to contradict 
previous studies showing no downcore variations in 
8 3 C of particulate organic carbon in sediments from 
Cape Lookout Bight [Blair et al., 1987; Martens et 
al., 1992; Blair and Carter, 1992]. However, the 
quantity of organic carbon converted to C0 2 and 
methane during the 1 14-day experiment (40.4 mM) 
represents a small fraction (~2%) of the particulate 
organic carbon pool. Therefore it is not surprising 
that preferential remineralization of a 13 C-rich labile 
fraction is not detected in S l3 C profiles of total 
organic carbon. 

Since C0 2 and methane are both end products of 
organic matter decomposition in sulfate-depleted sed- 
iments, the isotopic composition of organic carbon 
remineralized during the period of methane produc- 
tion is calculated as the weighted average 8 13 C of the 
ZC0 2 and methane produced: 


S 13 C-Cor G [MP] = 

(5 13 C-ZCQ 2 ) FIN [ZCQ 2 ] fin 

[IC0 2 ] F in - [£C0 2 ]i NI + [CH 4 ] fin - [CH 4 ] ini 

(8 C-ZC0 2 ) ]NI [£CQ 2 ]ini 

[ZC0 2 ] FIN - [ZC0 2 ] INI + [CH 4 ]p| N - [CH 4 ]ini 

(8 13 C-CH 4 )p] N [CH 4 ] fin 

[XC0 2 ]p IN - [IC0 2 ] rN i + [CH 4 ] F1N - [CH 4 ] ini 

(S i3 C-CH 4 ) IN i [CH 4 ] IN i 

[EC0 2 ]pin - [ZC0 2 ]ini + [CH 4 ] FIN - [CH 4 ] jjsjj 

where INI and FIN denote the concentration or 8 13 C 
value at the start of the methane production period 
(day 42) and at the end of the experiment (day 1 14), 
respectively. The ZC0 2 and methane concentrations 
at day 42 were 55.2 and 0. 12 mM, respectively 
(Figure 2). The final methane concentration 
(13.5 mM, Figure 6) was calculated by correcting for 
gas loss using (1). The final XC02 concentration 
(64.6 mM) was also corrected for gas loss using an 
equation analogous to (1). The ZC0 2 production rate 
was set equal to 0.7 times the total methane pro- 
duction rate (MPR). This stoichiometry is based on 
the oxidation state of labile organic carbon inferred 
from the AXC0 2 /AS0 4 2 ratio (-1.7). The 8 13 C-CH 4 
and 8 13 C-IC0 2 values at day 42 were -61.4 %o and 
-14.1 %o, respectively (Figure 5). The final 8 13 C- 
CH 4 value (-66.6 %o) was corrected for gas loss by 
integrating 8 13 C-MP (equation (5)) over the entire 
1 14-day interval. Note that this value represents the 
cumulative 8 l3 C-CH 4 and differs substantially from 
the measured value at day 1 14 (-50.8 %o). The final 
8 13 C-ZC0 2 value (-3.9 %o) was also corrected for 
gas loss using an equation analogous to (5). How- 
ever, for the ZC0 2 pool, isotopic fractionation 
between ZC0 2 (aq) and C0 2 (g) [Deuser and Degens, 
1967; Wendt, 1968] was taken into account. The 
8 C value of organic carbon remineralized during 
the period of methane production (S 13 C-C ORG [MP]) is 
-16.0 %o (equation (7)). 

The isotope mass balance provides a check on the 
accuracy of the 8 13 C-MP calculation. The 5 I3 C value 
of organic carbon remineralized during the period of 
sulfate reduction (8 13 C-Corg[SR]) is based entirely 
on measured concentration and 8 13 C values. In 
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contrast, the 8 13 C value of organic carbon reminer- 
alized during the period of methane production 
(5 13 C-Corg[MP]) is based on concentration and 8 13 C 
values corrected for gas loss according to (2)-(5). 
Agreement between 8 13 C-Corg[SR] (-15.6 %o) and 
5 i 3 C-C 0 rg[MP] (-16.0 %o) indicates that calculated 
5 13 C-MP values are consistent with the constraint 
that the microcosm maintain an overall isotope mass 
balance. 

8 1< C-CH 4 From C0 2 and Acetate 

The relationship between 8 13 C-CH 4 and the 
methanogenic pathway can be complicated by con- 
current changes in the isotopic composition of the 
methanogenic substrates (C0 2 and acetate) and/or 
variability in the magnitude of the kinetic isotope 
effect for the methane production pathways (CO2 
reduction and acetate fermentation). In this section, 
we examine how these factors may have played a 
role in controlling the stable carbon isotopic compo- 
sition of methane produced by the sediment 
microcosm. 

Methane from C0 2 reduction. There is a large 
kinetic isotope effect associated with bacterial C0 2 
reduction to methane. The magnitude of this effect is 
expressed as an isotope fractionation factor (occr), 
defined as the ratio of methane production rates from 
12 C- and 13 C-C02 normalized to the isotopic abun- 
dance in the reactant molecule [Rees, 1973]: 


_ I2 C-MPC/[‘ 2 CQ 2 ] 
CtcR_ 13 C-MPC/[ 13 C0 2 ] ’ 


( 8 ) 


where MPC denotes the rate of methane production 
via C() 2 reduction. 

If we let R M pc = 13 C-MPC/ 12 C-MPC and 
Ricoj] = [ 13 C0 2 ]/[ 12 C0 2 ], (8) can be written as 


Rmpc=^«-. W 

It follows from the definition of 8 I3 C [Craig, 1957] 
that 


8 13 C-MPC = 


8 13 C-C0 2 + 10 3 

OCR 


10 3 , 


(10) 


where 8 13 C-MPC represents the isotopic composition 
of methane produced from CO2 reduction. This treat- 
ment assumes that COo (as opposed to HCO3 or 


C0 3 2 ') is the form of inorganic carbon utilized by 
methanogenic bacteria [Fuchs et al, 1979]. 

The isotopic composition of methane produced 
from C0 2 at any time during the magnum syringe 
experiment can be determined from the 8 l3 C of the 
C0 2 pool and the isotope fraction factor (equation 
10). While 8 13 C-C0 2 values are easily calculated 
from 8 13 C-IC0 2 data [Deuser and Degens, 1967], 
there is considerable uncertainty regarding the iso- 
tope fractionation factor for methane production from 
C0 2 . Laboratory cultures of methanogenic bacteria 
grown autotrophically on H 2 -C0 2 yield fractionation 
factors that range from 1.025 to 1.061 (mean: 1.041, 
n = 7) [Oremland, 1988]. It is not clear whether this 
high degree of variability is due to differences in 
experimental conditions or differences between 
methanogenic species. For example, a single species 
(Methanobacterium thermoautotrophicum) grown in 
a sealed batch culture [Games and Hayes, 1976] and 
a flow-thru fermentation system [Fuchs et al., 1979] 
yielded isotope fractionation factors that differed by 
10 °/oo. Likewise, different species ( Methanobac- 
terium barkeri and Methanosarcina bryantii) grown 
under identical conditions had isotope fractionation 
factors that differed by 28 %o [Games and Hayes, 
1976], 

Aside from interspecies and intraspecies vari- 
ability, there are reasons why in vitro fractionation 
factors may not apply to natural methanogen popu- 
lations. First, the magnitude of the fractionation 
factor may be a function of the methane production 
rate. Laboratory cultures fed an 80/20 ( v /v) mixture of 
H 2 -C0 2 generally have higher production rates than 
methanogenic bacteria in natural ecosystems. 

Second, most laboratory measurements of isotope 
fractionation by methanogenic bacteria were con- 
ducted at elevated temperatures (37° to 65°C). 
Thermodynamics predicts that carbon isotope 
fractionation between C0 2 and methane decreases 
with increasing temperature [Bottinga, 1969], 
although it is not clear whether the kinetic isotope 
effect for bacterial methane production follows the 
same relationship. Third, the isotope fractionation 
factor may be species-specific and the bacteria used 
in the culture studies are unlikely to be representative 
of natural methanogen populations. 

One approach to deriving fractionation factors 
applicable to methanogenic populations in aquatic 
sediments is based on in situ isotopic differences 
between methane and C0 2 . In sediment where C0 2 
reduction is the exclusive methane source, the 
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isotope fractionation factor can be estimated from 
paired 5 13 C-C02 and 5 I3 C-CH 4 data (equation (10)). 
Evidence based on the deuterium content of methane 
and its formation water suggests that C0 2 reduction 
is the dominant methanogenic pathway in many 
marine sediments [Whiticar et al., 1986]. An 
extensive compilation of paired 5 13 C-CH4 and 
8 l3 C-C0 2 data from a wide variety of marine sedi- 
ment environments predicts a fractionation factor for 
C0 2 reduction of 1.070 ±0.020 [Whiticar et al., 
1986]. 

A more rigorous approach to estimating 
fractionation factors for C0 2 reduction in aquatic 
sediments takes into account the contribution of 
different methane production pathways, the intra- 
molecular isotopic composition of acetate, and the 
kinetic isotope effect for acetate fermentation. This 
approach was used to estimate fractionation factors 
for C02 reduction in Cape Lookout Bight sediments 
[Blair etal., 1992]. Isotope fractionation factors for 
summer (1.056 ±0.004) and winter (1.062 ±0.002) 
conditions fall within the range of values calculated 
from the paired 8 13 C-CH 4 and 8 13 C-C0 2 data. 

The isotope fractionation factors for Cape 
Lookout Bight sediment [Blair et al., 1992] were 
used to calculate S 13 C-MPC values (equation (10)) 
for methane produced by the sediment microcosm. 
An average Ocr value (1.06) was used because 
seasonal variations are relatively small. A rather large 
range (±0.01) was adopted to allow for uncertainty 
in the value of the isotope fractionation factor. 

The stable carbon isotopic composition of 
methane derived from C0 2 is denoted by the stippled 
region in Figure 7; the upper and lower limits corre- 
spond to Ocr = L05 and Ocr = L07, respectively. 
The trend toward more positive 8 l3 C-MPC values 
with time is due to the shift in S ]3 C-£C0 2 during the 
period of methane production (Figure 5a). The pro- 
gressive enrichment of 13 C in the XC0 2 pool is due, 
at least in part, to the kinetic isotope effect associated 
with C0 2 reduction: l2 C0 2 is converted to methane at 
a faster rate than 13 C0 2 , leaving the residual XC0 2 
pool isotopically heavier. 

The methane produced between 50 and 75 days 
is depleted in 13 C compared to the methane produced 
by C0 2 reduction. The methane produced after 
85 days is enriched in 13 C relative to the methane 
derived from C0 2 . Although we cannot rule out the 
possibility that the isotope fractionation factor for 
C0 2 reduction varied over the course of the experi- 
ment, CXcr values substantially less than 1.05 or 


greater than 1.07 are unlikely. Production of methane 
with S 13 C values well outside the range expected for 
CCb reduction are in all likelihood related to methane 
production from acetate. 

Methane from acetate fermentation. The S 13 C of 
methane produced from acetate (8 l3 C-MPA) cannot 
be calculated directly (using an equation analogous to 
equation (10)) because the 8 13 C of acetate in the 
magnum syringe sediment is unknown. However, 
8 l3 C-MPA can be estimated by an isotope mass 
balance calculation: 

5»c-mpa. (8 ' !c - m1>) -^- F)5 ' 1c - mpc , (11, 

where F denotes the fraction of methane derived 
from acetate (Figure 4b). The large range in 
8 13 C-MPC values, due to uncertainty in the isotope 
fractionation factor for C0 2 reduction, limits the 
accuracy of the S l3 C-MPA calculation. It is possible, 
however, to place constraints on the isotopic compo- 
sition of the methane produced by acetate fermen- 
tation. In the following section, we estimate 
8 13 C-MPA for two periods during the experiment: 

(1) the 51 to 58 day interval corresponding to mini- 
mum S l3 C-MP values and (2) the 86 to 93 day inter- 
val corresponding to maximum S 13 C-MP values. 

During the 51 to 58 day interval, the methane 
produced by the sediment microcosm was depleted in 
13 C (8 13 C-MP = -94 %o) compared to the methane 
produced from C0 2 (8 l3 C-MPC > -83 %») (Figure 
7). This implies a source of 13 C-depleted methane 
derived from a pathway other than C0 2 reduction. 
Acetate fermentation is the most likely source of this 
methane; the lack of methane production in the pre- 
sence of sulfate argues against significant methane 
production from noncompetitive substrates. Acetate 
fermentation accounted for 34% of the methane pro- 
duced at 58 days (Figure 4b). Isotope mass balance 
(equation (11)) requires that the methane produced 
by acetate fermentation during the 51 to 58 day 
interval have a 8 l3 C value less than -115 %o. This 
extreme level of 13 C depletion is only possible if the 
methyl group of the acetate molecule is highly 
depleted. 

Net acetate accumulation during the 5 1 to 58 day 
interval (Figure 2b) indicates that gross production 
exceeded consumption. Since acetoclastic methan- 
ogens were not substrate limited, isotopic fraction- 
ation of the acetate methyl carbon was not influenced 
by closed system effects. Under these conditions, the 
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carbon isotopic composition of the acetate methyl 
group (8 l3 C- + CH 3 COOH) can be calculated by 
rearranging an equation analogous to (10): 

8 13 C-*CH 3 COOH = 

S 13 C-MPAa AF + (a AF - 1) 10 3 , (12) 

where a A F is the isotope fractionation factor for 
acetate dissimilation. Laboratory estimates of 
range from 1.02 to 1.03 for pure cultures of the 
acetate-utilizing methanogen Methanosarcina barkeri 
[Risatti and Hayes, 1983; Krzycki et al., 1987]. 
Comparable values (1.032 ± 0.014) were estimated 
for methanogenic sediments from Cape Lookout 
Bight [Blair and Carter, 1992]. The latter values, 
inferred from measured downcore shifts in the 
isotopic composition of the acetate methyl carbon, 
are adopted for the following calculations. Given a 
8 I3 C-MPA of -1 15 %o for the 51 to 58 day interval 
(from equation (11)), we calculate a S I3 C of the 
acetate methyl carbon of -74 to -99 %» (equation 
(12)). The large range in the 8 l3 C-*CH 3 COOH 
values stems from uncertainty in a^. 

Autotrophic acetogenic bacteria provide a mecha- 
nism for producing acetate with a methyl group that 
is highly depleted in l3 C. Gelwicks et al. [1989] have 
shown that pure cultures of Acetobacterium woodii 
grown autotrophically on C0 2 and H 2 produce iso- 
topically homogeneous acetate that is depleted in l3 C 
by 59 %o relative to the IC0 2 pool. If the fraction- 
ation factor for Acetobacterium woodii in pure 
culture is applicable to microbial populations in the 
sediment microcosm, autotrophic acetogenesis 
during the 51 to 58 day interval would produce 
acetate with a methyl 8 13 C value of -7 1 %o. 

Although this 8 13 C-*CH 3 COOH value is slightly 
outside the range predicted above, it demonstrates the 
plausibility of extreme 13 C depletion in methane 
derived from autotrophically produced acetate. 

At 86 days, acetate fermentation accounted for 
87% of the methane produced by the sediment 
microcosm (Figure 4b). The methane produced 
during the 86 to 93 day interval was enriched in 13 C 
(S 13 C-MP = -42 °/oo) compared to the methane 
derived from CO 2 reduction (8 I3 C-MPC< -62 %o) 
(Figure 7). Isotope mass balance constraints indicate 
that methane produced by acetate fermentation during 
this period had a 8 13 C value of -37 to -39 %x> 
(equation (11)). 

The acetate pool was nearly quantitatively con- 
sumed during the 86 to 93 day interval (Figure 2a). 


Assuming no isotopic fractionation during methyl 
group oxidation [Blair and Carter, 1992], the rela- 
tionship between 8 13 C-*CH 3 COOH and 8 13 C-MPA 
is given by a closed- system Raleigh equation 
[Hayes, 1983]: 

8 13 C-*CH 3 COOH = 

8 13 C-MPA - ((Xa - [In (1-/)] 10 3 , (13) 

a A F / 

where / is the fraction of the acetate methyl carbon 
pool converted to methane. The / value (0.80) was 
calculated from the rate of methane production from 
acetate at 86 days (Figure 3b) divided by the net 
acetate depletion rate during the 86 to 93 day interval 
(Figure 2b). From (13), we calculate 8 13 C- 
*CH 3 COOH values ranging from -24 to -33 %o. 
Although the range of calculated values is large 
(reflecting uncertainty in a^), they overlap with 
values expected for acetate produced by fermentative 
processes (-21 to -25 %o) [Blair and Carter, 1992]. 

The 8 l3 C of methane produced from acetate 
(calculated by isotope mass balance, equation (11)) 
varied from less than - 1 15 °/oo (51 to 58 day interval) 
to greater than -39 %o (86 to 93 day interval). Such 
large variations cannot be accounted for by changes 
in the isotope fractionation factor for acetate dissim- 
ilation. Rather, the variations in 8 l3 C-MPA imply 
changes in the isotopic composition of the acetate 
methyl carbon. 

Microbial Controls on 8 }3 C-CH 4 

We have shown that large variations in the 8 l3 C 
of methane produced by the sediment microcosm 
were driven by shifts in the relative rates of C0 2 
reduction and acetate fermentation, as well as 
changes in the isotopic composition of the C0 2 pool 
and the acetate methyl carbon. In this section, we 
present a scenario describing the microbial processes 
that may have played a role in controlling the vari- 
ations in isotopic composition of methane produced 
by the microcosm. 

Organic matter degradation in the sediment 
microcosm can be divided into three distinct periods: 
(1) the sulfate reduction period, (2) the sulfate- 
methane transition period, and (3) the methane 
production period. During the sulfate reduction 
period (0 to 42 days), high sulfate reduction rates 
(Figure 3a) were supported by rapid degradation of 
relatively fresh organic matter. Sulfate concenn*ations 
(Figure 2a) throughout this period exceeded the level 
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at which sulfate-reducing bacteria become sulfate- 
limited (0.2 mM) [Ingvorsen et al., 1984], In the 
presence of nonlimiting sulfate concentrations, 
sulfate-reducing bacteria maintain acetate and H 2 
concentrations at a level too low to support growth of 
methanogenic bacteria [Lovley et al., 1982; Lovley 
and Klug, 1986], The acetate and H 2 produced by 
organic matter degradation were essentially unavail 
able to methanogenic bacteria and total methane 
production rates were negligible (Figure 4a). 

During the sulfate-methane transition period 
(42 to 60 days), sulfate concentrations dropped to 
threshold levels for uptake by sulfate-reducing 
bacteria (10 to 20 pM) [Ingvorsen et al., 1984], 
Methane began to accumulate in the magnum syringe 
(Figure 2a), suggesting that methanogens had 
replaced sulfate reducers as the primary acetate and 
H 2 consumers. Methane production was initially 
dominated by C0 2 reduction (Figure 3b), possibly 
because methanogens growing on H 2 -C0 2 have 
faster doubling times than acetoclastic species 
[Vogels et al., 1988]. The first methane produced by 
the microcosm had a 8 13 C of -80 °/oo, within the 
range of values expected for methane produced from 
CO, (Figure 7). 

The accumulation of acetate during the sullate- 
methane transition period (Figure 2b) suggests that 
the existing population of acetate-utilizing methano- 
genic bacteria was unable to keep pace with acetate 
production. It is not known whether the H 2 reservoir 
experienced a similar accumulation, but relatively 
slow rates of C0 2 reduction during the sulfate- 
methane transition (Figure 3b) suggest that the 
population of methanogenic bacteria was not fully 
developed. These conditions would be favorable for 
H 2 -consuming acetogenic bacteria that normally 
cannot compete with sulfate-reducing or methano- 
genic bacteria [Lovley and Klug, 1983]. Most 
acetogenic bacteria capable of growing on H 2 -C() 2 
are also able to grow by fermenting sugars to acetate 
[Dolfing, 1988]. Therefore acetogenic bacteria 
(unlike methanogens) can exist in a commensal 
relationship with sulfate-reducing bacteria, and 
would be well poised to take advantage of a tem- 
porary reduction in H 2 demand that might follow 
sulfate depletion. Acetate produced from H 2 -CO, by 
autotrophic acetogenic bacteria would be highly 
depleted in I3 C [Gelwicks et al., 1989], Subsequent 
methane production from this acetate could account 
for the l3 C-depleted methane (-94 °/oo) produced at 
the end of the transition period (Figure 7). 


The methane production period (60 to 1 14 days) 
can be divided into two subperiods. During the first 
subperiod (60 to 90 days), methane production was 
dominated by acetate fermentation (Figure 3b). The 
high methane production rates were fueled by the 
large acetate pool that had begun accumulating during 
the sulfate-methane transition (Figure 2b). Methane 
production rates from C0 2 reached a plateau during 
this subperiod (Figure 3b), suggesting that H 2 - 
consuming methanogens were keeping pace with Re- 
production. Since acetogens cannot successfully 
compete for H 2 with a well-developed methanogen 
population [Lovley and Klug, 1983], the acetate that 
accumulated throughout most of this subperiod was 
probably synthesized via the sugar fermentation path- 
way. The methyl carbon in acetate produced by sugar 
fermentation would be isotopically "heavy," and thus 
give rise to methane enriched in l3 C [Whiticar et al., 
1986]. The l3 C-rich methane (-42 %o) produced 
toward the end of the initial methane production 
subperiod (Figure 7) is consistent with the predom- 
inance of methane production from acetate synthe- 
sized by the fermentative pathway. 

The second methane production subperiod (90 to 
1 14 days) marks the return to a stable microcosm. 
Acetate concentrations were depleted to steady state 
levels (Figure 2b) and total methane production rates 
decreased gradually (Figure 4a) as the abundance of 
labile organic matter declined. Assuming that 
methanogens were the exclusive acetate and H 2 sinks 
during this period, the relative availability of acetate 
and H 2 must have been controlled by their respective 
production rates. Methane production was dominated 
by C0 2 reduction (Figure 3b), suggesting that H 2 
was produced in stoichiometric excess of acetate. 

The methane produced during the final subpenod 
was enriched in 13 C (-46 to -60 %o) relative to the 
methane produced during the sulfate-methane tran- 
sition (-80 to -94 %o). Although CO, reduction 
accounted for -70% of the methane produced during 
both these periods (Figure 4b), the methane pre- 
cursors were enriched in 13 C during the latter period. 
The ICCh pool was progressively enriched in C 
following sulfate depletion (Figure 5a), presumably 
due to preferential uptake of 12 C() 2 by methanogenic 
and possibly autotrophic acetogenic bacteria. The C 
enrichment in the acetate methyl group may be related 
to changes in the pathway of acetate synthesis 
(i.e., autotrophic acetogenesis during the sulfate- 
methane transition; sugar fermentation during the 
methane production period). 
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Comparison of Laboratory Microcosm and Natural 
Sediment 

A laboratory microcosm is a highly artificial 
system. Removing sediment from its natural setting 
and placing it in a sealed container will inevitably 
lead to perturbations. The results of the microcosm 
experiment must be compared with in situ data to 
evaluate whether the processes that control isotopic 
variability in the sediment microcosm may also 
operate in undisturbed marine sediment. 

The stable carbon isotopic composition of 
methane contained in gas bubbles emitted from Cape 
Lookout Bight sediments (Figure 8 a) varies by 
~10 %o over the course of a year [Martens et al., 
1986; Burke et al„ 1987], Between February and 
May, 8 13 C-CH 4 values remain relatively constant at 
-63 to -65 °/oo. For a short period at the end of May, 
^ C-CH 4 values drop to -68.5 °/oo. The transient 
decline is followed by a shift toward more positive 
5 1 C values that continues through the middle of 
August. Between August and November, 8 13 C-CH 4 
values gradually decline from -57 to -62 %o. 

Martens et al. [1986, p. 1300] hypothesize that the 
observed seasonal isotopic variations are ... 
controlled by changes in pathways of methane 
production and cycling of acetate and molecular 
hydrogen within the sediments". 

There is a dramatic increase in the acetate inven- 
tory in Cape Lookout Bight sediments between May 
and August (Figure 8 b). The bulk of the summertime 
acetate inventory is contained in a narrow subsurface 
peak centered at 10 to 15 cm [Sansone and Martens, 
1982], The vertical extent of the sulfate reduction 
zone in Cape Lookout Bight decreases during the 
spring, as sulfate reduction rates increase in response 
to higher sediment temperatures [Crill and Martens, 
1987], The appearance of the subsurface acetate peak 
at the end of May corresponds to the sulfate- 
depletion depth crossing the 15 cm horizon [Sansone 
and Martens, 1982], 

Sediment accumulation rates for Cape Lookout 
Bight (determined from 7 Be inventories) indicate that 
sediment residing at the 10 to 15 cm interval in May 
was deposited at the sediment surface during the late 
fall [Canuel et al., 1990], As a result of relatively 
cool sediment temperatures, organic matter deposited 
in the fall is effectively preserved through the winter 
and spring, reaching the 10 to 15 cm horizon with a 
large fraction of its labile component intact. As sedi- 
ments warm during the spring, this labile material is 
able to support high rates of acetate production. 
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Fig. 8 . Seasonal changes in Cape Lookout Bight 
sediment, (a) 8 13 C of methane in gas bubbles 
collected during 1983 to 1984 (open circles denote 
data from Martens et al. [1986]; closed circles denote 
data from Burke et al. [1987]); (b) acetate inventory 
calculated by integrating acetate concentration-depth 
distributions (0 to 36 cm) determined on sediment 
cores collected in 1990. The error bars in Figure 8 a 
denote the standard deviation based on two to four 
replicate samples. 


The processes that lead to rapid acetate accu- 
mulation in the sediment column (Figure 8 b) and the 
sediment microcosm (Figure 2 b) are presumably 
similar, high rates of acetate production accompanied 
by a sudden reduction in the rate of acetate con- 
sumption. In the case of the sediment microcosm, the 
decline in acetate consumption appears to be caused 
by the cessation of sulfate reduction (due to sulfate- 
limitation) and a lag in acetate consumption by 
methanogenic bacteria (due the slow growth rate of 
acetoclastic methanogens). 

There is a striking similarity in the relationship 
between 8 13 C-CH 4 and acetate concentration (or 
inventory) for both the microcosm and the undis- 
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turbed sediment (cf. Figures 2b, 7, and 8). Foi both 
systems, methane produced (or emitted) at the 
beginning of the acetate accumulation period is rela- 
tively depleted in l3 C. Likewise, methane produced 
(or emitted) during the period of rapid acetate 
depletion is relatively enriched in 13 C. Seasonal 
variations in 6 l3 C-CH 4 (~ 10 %») are much smaller 
than variations in the microcosm (-50 %o), pre- 
sumably due to inertia imposed by the sediment 
column methane inventory. 

The similarities between laboratory and field data 
suggest that similar processes control the isotopic 
composition of methane produced by the laboratory 
microcosm and natural sediment. The shoaling of the 
sulfate-depletion depth to the sediment horizon con- 
taining labile organic material may initiate a scenario 
similar to that described for the sediment microcosm. 
First, sulfate depletion may lead to a temporary 
reduction in competition for acetate and H 2 , resulting 
in acetate pooling and possibly, pioduction of C- 
depleted acetate by autotrophic acetogenic bacteria. 

The 13 C-depleted methane emitted at the end of May 
is consistent with a predominance of the C0 2 - 
reduction pathway (due to the slow growth rate of 
acetoclastic methanogens) and possibly, methane 
production from 13 C-depleted acetate. Second, 
acetate fermentation may become the dominant 
methane production pathway as the acetoclastic 
methanogens respond to the elevated acetate concen- 
trations. Since acetogens cannot compete for H 2 with 
an active methanogen population, autotrophic aceto- 
genesis is unlikely. The 13 C-rich methane emitted 
from June thru August is consistent with acetate 
dissimilation being the dominant methane production 
pathway. Third, the methane production pathway 
may shift in favor of C0 2 reduction after the acetate 
inventory is fully depleted. The shilt toward C- 
depleted methane from September thru November is 
consistent with increased importance of C0 2 
reduction after depletion of the acetate reservoir. 

CONCLUSIONS 

The laboratory microcosm experiment demon- 
strates that a single sample of anoxic marine sediment 
incubated at constant temperature can produce 
methane with 8 13 C values that vary from -94 to 
-42 %o. The variations in 8 13 C-CH 4 were caused by 
changes in the relative importance of acetate fer- 
mentation and Ct)i reduction, as well as changes in 
the isotopic composition of the methane precursors 
(acetate and C0 2 ). Concurrent methane oxidation and 


production were observed, but oxidation was a small 
fraction (2.4%) of the gross production. It is unlikely 
that methane oxidation was an important factor in 
controlling the isotopic composition of methane. 

Variations in the relative importance of C0 2 re- 
duction and acetate fermentation can be understood in 
terms of changes in microbial processes that occurred 
as sediment bacteria responded to sulfate depletion. 

The initial predominance of the C0 2 reduction path- 
way is consistent with the relatively slow growth rate 
of acetate-utilizing methanogens. The subsequent 
shift in pathway to acetate fermentation appears to be 
a response to the elevated acetate concentrations that 
had accumulated following sulfate depletion. The 
return to C0 2 reduction as the dominant pathway 
toward the end of the microcosm experiment corre- 
sponds with the rapid and nearly quantitative deple- 
tion of the acetate pool. 

The apparent changes in the isotopic composition 
of the acetate methyl carbon may be related to 
changes in the pathway of acetate synthesis, which in 
turn, were controlled by competition for H 2 . 

Changes in isotopic composition of the Z.C 0 2 pool 
during the microcosm experiment are presumably 
due to isotopic fractionation associated with C0 2 
reduction by methanogenic and possibly autotrophic 
acetogenic bacteria. 

A comparison of microcosm and field data 
reveals a striking similarity in the relationship 
between 8 l3 C-CH 4 and acetate concentration or 
inventory. This suggests that the processes that 
control in situ seasonal variations in the isotopic 
composition of methane emitted from Cape Lookout 
Bight sediments may be similar to those that con- 
trolled variations in the isotopic composition of 
methane produced by the sediment microcosm. This 
experiment demonstrates that laboratory-based 
process studies coupled with field measurements 
provide a powerful tool for interpieting and 
understanding stable isotope signatures in nature. 
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